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ABSTRACT 

The centers of stellar spheroids less luminous than ~ 10 10 L© are often marked by the presence of 
nucleated central regions, called "nuclear star clusters" (NSCs). The origin of NSCs is still unclear. 
Here we investigate the possibility that NSCs originate from the migration and merger of stellar 
clusters at the center of galaxies where a massive black hole (MBH) may sit. We show that the 
observed scaling relation between NSC masses and the velocity dispersion of their host spheroids 
cannot be reconciled with a purely "in-situ" dissipative formation scenario. On the other hand, the 
observed relation appears to be in agreement with the predictions of the cluster merger model. A 
dissipationless formation mode l also reproduces the observed relation between the size of NSCs and 
their total luminosity, R oc v / £nsc- When a MBH is included at the center of the galaxy, such 
dependence becomes substantially weaker than the observed correlation, since the size of the NSC 
is mainly determined by the fixed tidal held of the MBH. We evolve through dynamical friction a 
population of stellar clusters in a model of a galactic bulge taking into account dynamical dissolution 
due to two body relaxation, starting from a power-law cluster initial mass function (CIMF) and 
adopting an initial total mass in stellar clusters consistent with the present-day cluster formation 
efficiency of the Milky Way (MW). The most massive clusters reach the center of the galaxy and merge 
to form a compact nucleus; after 10 10 years, the resulting NSC has properties that are consistent with 
the observed distribution of stars in the MW NSC. When a MBH is included at the center of a galaxy, 
globular clusters are tidally disrupted during inspiral, resulting in NSCs with lower densities than 
those of NSCs forming in galaxies with no MBHs. We suggest this as a possible explanation for the 
lack of NSCs in galaxies containing MBHs more massive than ~ 10 8 M© . Finally, we investigate the 
orbital evolution of globular clusters in giant elliptical galaxies which are believed to always host a 
MBH at their center rather than a NSC. In these systems an additional mechanism can prevent a 
NSC from forming: the time for globular clusters to reach the center of the galaxy is much longer 
than the Hubble time. 

Subject headings: galaxies: Milky Way Galaxy- Nuclear Clusters - stellar dynamics 



1. INTRODUCTION 

Many galaxies, over the whole Hubble sequence, 
show nucleated central regions often referred to 
as "nuclear star clusters" (NSCs). These systems 
are among the densest star clusters observed, with 
effective radii of a few pars ecs and central lumi- 
nosities up to ~ 10 7 L P) (Mat thews fe Gallagher I 
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There is no consensus on how NSCs form. A dissipa- 
tive origin bases on the hypothesis of radial gas inflow 
into the galactic center and requires efficient d issipation 
mechanisms to work (e.g., lLoose et al. II1982T) . In this 
model a NSC consists most l y of stars that formed lo- 
cally (ISchmnerer et al. 1 120061 [20081: IMilosavlievic 112004: 
Emsellem fc van de Ven I 120081 : Shlosman & Begclman 
1989UBekkill2007D . 

Alternatively, NSCs could have a dissipationless ori- 
gin in which massive stellar (globular-like) clusters mi- 
grate to the center due to dynamical friction and merge 
to form a den se nucleus growing up to the size of ob- 
served NSCs (jTremaine et al. I fl975). Observations of 
NSCs in dE galaxies suggest that the majority of these 
systems could be the result of accumulating mass in the 
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form of orbitally decayed globular clusters. Numerical 
studies have also shown that such a formation model is 
consist ent with the measured sizes and luminosities of 
nuclei dBekki et al71l2004J: iCapuzzo-Dolcetta fe Miocchi I 
120081 : lHartmann et al. 1 1201 lh . In addition, stellar popu- 
lation synthesis studies of NSC spectra suggest that most 
of the mass usually resi des in stars as old as typical glob - 
ular clusters (~ 10 Gvr: iFiger et"ai~|[200l lB5ker If^OlOr) . 
On the other hand, the observed correlation between col- 
ors and luminosities together with the complex star for- 
mation histories that often characterize the central re- 
gion of galaxies, may be difficult to explain on the basis 
of a dissipationless origin, unless there is some contribu- 
tion from continuous or recurrent star fo rmation in addi- 
tion t o the ancient globular cluster stars ([Antonini et al. I 
|20H . 

Due to their small sizes and their crowded stellar fields, 
galactic nuclei beyond the Local Group are typically un- 
resolved, and the only quantities that can be determined 
are integrated properties such as half-mass radius and 
total luminosity. A radial density profile and velocity 
structure can be rel iably determined on l y for the Milky 
Way (MW) NSC (IGenzel et al. I [20(51 I Schodel et al. I 
[20071: [Graham fc Spitler 1 120091: IQh et al. N2009IL which, 
due to its pro ximity (~ 8 kpc) , can be resolved into ind i- 
vidual stars (|Schodel et al. 1 120071: ISchodel et al. 1120091) . 
In addition to the MW, NGC 205 also has a spatially- 
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resolved NSC (Figure 2 of lMerritt l2Q p) . The MW NSC 
has an estimated mass of ~ 10 7 M Q (jLaunhardt et al."1 
12001 iSchodel et aD l2008h and hosts a massive black 
hole (MBH) who se mass, ~ 4 . 3 x 10 6 Mp), is uniquely 
well determined ([Genzel et al. 1 120031 : iGhez et al. 1120081 : 
IGillessenl I2009D . A number of ot her galaxie s also 
contain both a NSC and a MBH (iSeth et aTl 120081: 
iGraham fe Driver ll2007tlGraham fe Spitler 1120091 ). that 
have comparable masses. In models of NSCs, the dy- 
namical influence of a MBH should therefore be consid- 
ered, at least in bulges brighter than a bout 10 9 L ( r l which 
are b elieved to always contain a MBH ([Ferrarese fe Fordl 
120051) . 

More recently, we heave presented large-scale A-body 
simulations of the inspiral and merger of massive clus- 
ters i n the inner regions of the Galaxy (jAntonini et al. 1 
I2012D . We showed that current observational constraints 
are consistent with the hypothesis that a large fraction 
of the MW NSC mass is in old stars brought in by in- 
falling globular clusters. In this paper, we expand on 
this previous work and use simple analytical considera- 
tions to investigate the possibility that globular clusters 
can migrate through dynamical friction in the center of 
galaxies and form compact stellar nuclei. In particular, 
we focus on how the presence of central MBHs at the 
center of galaxies can impact the merger hypothesis for 
the formation of their NSCs. 

This paper is organized as follows. We start in §2 
by comparing some of the observed scaling relations for 
NSCs with the same relations predicted in both the 
merger and the gas model. In §3 and §4 we develop a 
simple analytical model for studying the orbital decay of 
globular clusters in galaxies like the MW, and explore the 
effect that a central MBH in the galaxy has on the prop- 
erties of the resulting NSC. In §5 we use an approach 
similar to that of §3 to follow the orbital evolution of 
globular clusters near massive MBHs in the central re- 
gions of bright elliptical galaxies. In §6 we discuss a 
variety of astrophysical implications of a dissipationless 
origin of NSCs and conclude in §7. 

2. SCALING RELATIONS 

The study of NSCs has revealed a number of corre- 
lations between their masses and several global prop- 
erties of their host galaxies, such as velocity disper- 
sion and bulge mass. The existence of such correla- 
tions might indicate a direct link among large galac- 
tic spacial scales and the much smaller scale of the 
nuclear environment. While it was for long believed 
that NSCs and MBHs followe d similar scaling rela- 
tions with their host galaxies (|Ferrarese et al. I 120061 : 
IWehner fc Harris! [2006) , it is now well established that 
NCs do n ot follow any of the scaling relations defined 
by MBHs ([Balcells et al. If2007i: IGraham fc Spitler 112009b 
Graham I \2012\ iScott fc Graham I 120121 iLeigh et al. I 



2012ft . Observations also reveal that, unlike globular 



clusters, NSC half-light radii are lumi nosity dependent , 
increasing with inc reasing total mass ([Cote et al. ^ 2006: 
IForbes et al. 1 12008). This relation might contain impor- 
tant information on the processes that shaped the central 
regions of galaxies and their NSCs. 

We start here by comparing such observed correlations 
with predictions from both the dissipative and dissipa- 



tionless formation models. 

2.1. Mnsc — cr relation 

Of the global-to-nucleus relations, the most frequently 
referred to is the tight correlation between NSC mass, 
MmsCi and the host galaxy's velocity dispersion, a. 
iFerrarese et al. I (1200 6 ■ argued that such a correlation has 
a slope which is consistent with that of the well-known 
M. — a relation obeyed by the MBH mass, M # . More 
recently, a series of papers reached a different conclusion, 
suggesting that the NSC scaling relations are instead sub- 
stantially shal lower than the corresponding MBH scal- 
ing relations (IGraham fc Spitler 1 120091: IGraham 1 120121 : 
IScott fc Graham ll2012HLeigh et al. II2012D The v ersion 
of the Mnsc — °~ relation given in " IGraham I (|20T1) is 



log (6.83 ±0.07) 



(1) 



+ (1.57 ± 0.24) log(cr/70km s" 1 ) 



We mention here that the Mnsc — G relation might 
be a non-primary correlation, instead resulting from a 
projection of the fundamental plane, given the observed 
correlation between Mnsc and the total host galaxy- 
mass (M gx ): Mnsc oc M g 1 x 18±016 (e.g., see ILeigh et al I 

nnn). 

2.1.1. Dissipationless model 

A predicted Mnsc — °~ relation can be easily derived in 
the globular cluster merger model if we assume that the 
globular cluster distribution initially follows the stellar 
light and by using an isothermal sphere density model: 
p(r) — a 2 /2irGr 2 , where a is the ID velocity dispersion 
and G the gravitational constant. The cumulative mass 
within r is M(r) = rv 2 /G = 2ra 2 /G, with v c the velocity 
of a circular orbit. The dynamical friction coefficient for 
a stellar cluster of mass m r ] moving in an isothermal 
sphere model is ([Chandrasekhar 1119431 ): 
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with In A the Coulomb logarithm and X — v/y/2<r. Not- 
ing that 
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where L(= \J GM(r)r) is the orbital angular momentum, 
for a cluster moving on a circular orbit we obtain 
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Integrating Equation ([4]) yields 
i' in = f v^GrndlnA 



erf m _ ^L-* 2 
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Clusters initially within r; m reach the center at a time 
< t. 

Assuming that the total mass accumulated in the cen- 
ter is equal to the total mass in globular clusters that are 
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initia lly within rj n , then we obtain (e.g.. lTremaine et al. I 
1975) ' 



Af NSC = 2 5/4 (hi A 
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(6) 



where m* is the mass of the field stars, (m c \) is the av- 
erage globular cluster mass, and / is the initial number 
fraction of globular clusters (after galaxy formation) to 
the total number of stars in the galaxy. Since the term in 
square brackets in Equation (JB]) is a constant, the globu- 
lar cluster merger model predicts A/nsc « o -3 ^ 2 - 

In the MW, the initial total number of globular clus- 
ters, iV c i(0), can be recovered based on their observed lu- 
minosity function and semi-analytical modeling of mass- 
loss due to stellar evolution and due to tidal inte r action 
with the Galactic environment. iBonatto fc Bica 1 ()2012f ) 
inferred N c i(0) ~ 10 4 , which, compared to the total stel- 
lar mass of the Bulge, - 10 10 M©, yields / « 1CT 6 . We 
can rewrite Equation ([6]) as 



Mnsc = 3 x 10 7 M G 
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Despite its simplicity, our model reproduces the observed 
-Mnsc — o relation both in slope and normalization. 

2.1.2. Dissipative model 

iMcLaughlin et al. I (|2006t ) proposed a NSC in-situ star 
formation model regulated by momentum feedback . This 
model is an extension of the argument proposed bv lKingl 
(2003) to explain the M, — a relation, and invokes the 
formation of a massive nuclear cluster due to gas inflow 
and accumulation at the center of the galaxy during the 
early phases of galaxy evolution. Stellar winds and su- 
pernovae from a young nuclear cluster with a standard 
IMF produce an outflow with momentum flux given by 



II w AL Edd /c : 



XA-kGM] 
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(8) 



where ^Edd is the Eddington luminosity calculated by 
the stellar mass, c is the speed of light, A ~ 0.05 
is related to the mass fraction of massive stars, and 
k = 0.398cm 2 g -1 is the electron scattering opacity. The 
outflow is initially momentum conserving and produces 
an outward force on the gas in the bulge, whose weight 
is W(r) = GM gas {r)M(r)/r 2 , with M gas (r) the enclosed 
gas mass and M(r) the total enclosed mass of the galaxy. 
For an isothermal potential one finds 



(9) 



wher e f a 



0.16 



the baryonic mass frac- 
tion (jSpergel et al. II2003D . Requiring that the momen- 
tum output produced by the nuclear cluster balances the 



weight of the gas leads to the relation 
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A similar scaling relation can be derived in 
protog alaxies with non- isothermal dark matter 
haloes (|McQuillin fc McLaughli"n"ll2012h . 

This model is very attractive because it contains no 
free parameters. In addition, Equation (flOt is in good 
agreement with the M nsc — o~ correlation reported by 
iFerrarese etaD (|2006f ): 



1ok('^ )=(«.»] i 0.11. 



(11) 
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but, a s previously mentioned, it is at odds with lGraham I 
(2012) who finds NSC scaling relations considerably shal- 
lower than the corresponding MBH scaling relations. 
The difference between these two studies was due to 
3/2 the proper excl usion of nuclear disks in the sample of 
I Graham 1 ([2012) and the larger sample of NSCs used 
in this latter work. For these reasons we consider the 
results of Graham more robust and conclude that the 
(7) McLaughlin et al. model provides a poor description of 
the observed Mnsc — o correlation. This suggests that 
momentum feedback may be not relevant, which would 
be expected if the NSCs originated elsewhere and were 
subsequently deposited into their host galaxy centers. 

2.2. R — Mnsc relation 

The size of galactic nuclei clearly correlates with 
their luminosity, in the sense that brighter NSCs 
have la rger effective rad ii. The relation is approxi- 
mately (|C6te et al. II2006D 



R oc a/^NSC ) 



(12) 



where R is the NSC effective radius (or half-mass ra- 
dius) and L its total luminosity. This relation is consis- 
tent with the extrapolation of the R e f / — C relation for 
elliptical galaxies. 

2.2.1. Dissipationless model 

In the merger model the radius of the nucleus increases 
with increasing total mass as globular clusters merge. 
The brighter, larger mass nuclei are therefore pr e dicted 
to be spatially very extended. lAntonini et al. I (2012) 
used simple energy arguments to derive the size-mass 
relation for NSCs in galaxies with no central MBH. For 
the sake of completeness we repeat this simple calculation 
in what fallows. After a merger the NSC energy, Ef — 
—GMj/2Rf, equals the energy of the nucleus before the 
merger, Ei, plus the energy brought in by the globular 
cluster: 



Ef — Ei + E a 



E h 



(13) 



with Ei, the internal binding energy of the cluster, and 
E a —Gm c \Mi/2Ri its orbital energy before the merger. 
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The equations above permit expressing the mass, en- 
ergy, and radius of the nucleus recursively as 

(14) 
(15) 
(16) 



M j+1 = (j + l)M u 
jE j+1 = (j + l)E j +jE 1 , 



i) 2 RJ+i 



■■ i (j + 1) Rj 



Ri\j = 1,2,3,. 



where the subscript 1 denotes the initial nucleus, and, 
by assumption, A/i ~ m. Equations (I14j|16[) imply 
R oc Af 05 at the time the mass of the nucleus is still com- 
parable to the mass of the infalling clusters. After many 
mergers the nucleus is much more massive than one glob- 
ular cluster and the relation steepens to R oc M . After 
25 mergers, Equations (|14H16|) imply R « 5i?i . The half- 
mass radii of globular clusters are ~ 3 pc ([Jordan et al. I 
2005), irrespective of their luminosity, so for a nucleus 
assembled from 25 mergers, R ~ 15 pc. Such a value 
is in reasonable agreement with the measured half-mass 
radii for the brightest nuclei. We conclude that a model 
that attributes the origin of NSCs to the mergers of glob- 
ular clusters at the centers of galaxies is consistent with 
the sizes and luminosities of the nuclei. 

In a number of galaxies NSCs are observed to co- 
exist with MBHs. An example of such systems is 
the MW for which the mass o f the NSC, M NS c ~ 
10 7 M & (ILaunhardt et al. Il200l TSchodel et al. 1120081) . 
is somewhat comparable to the mass of the central black 
hole M . ~ 4 X 10 6 M Q (|Ghez et al. I [200l iGillessen I 
2009) . A MBH will disrupt clusters that pass within the 
radius 
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where ok is the one-dimensional velocity dispersion of 
a globular cluster, tk its core radius, ctnsc i s t ne ve ~ 
locity dispersion in the NSC and n n fl = GAf./er NSC 
is the influence radius of the MBH. In the presence of 
a MBH, the dependence of the half-mass radius R on 
Mnsc is substantia lly weaker than the observed correla- 
tion, R ~ M N f c (jAntonini et al. II2012T ). due to the fact 
that the size of the NSC is determined by the fixed tidal 
field from the MBH. When a stellar cluster is disrupted, 
stars that were initially within the cluster core will re- 
distribute locally at a distance rdi sr from the MBH. The 
density profile of the NSC will therefore have a core of 
characteristic size ~ rdi sr . 

2.2.2. Dissipative model 

Although the gas model remains somewhat more qual- 
itative there are some indications that the observed 
R — Mnsc relation might be difficult to reconcile with 
a purely d issipa tive scenario. 

Bckki (2007) performed fully self-consistent chemo- 
dynamical simulations to investigate how NSCs can 
form through dissipative gas dynamics. He found that 
compact nuclei can be formed via dissipative, repeated 
merger of gaseous (or stellar) clumps that develop from 
nuclear gaseous spiral arms due to local gravitational 
instability. These computations showed that fainter 
NSCs are likely to have a more diffuse configuration due 
to more negative feedback from SNe II which in turn 
can prevent NSCs from forming in faint galaxies. The 
fact that NSC formation is more-strongly suppressed by 



stronger feedback effects in less-luminous galaxies would 
explain why brigh ter dwarf galaxies (dE ) are more likely 
to contain NSCs (|van den Bergh Ill986l ). But this would 
also predict that more massive NSCs are less extended 
than their lower mass counterparts, which is the opposite 
of the observed trend. 

Further investigation of the above inconsistency is 
needed in order to determine how the results depend on 
the resolution limit of the simulations and on the adopted 
initial conditions. 

3. FORMATION OF NSCS 
3.1. Phase-space constraints 

The above discussion gives some level of reliability to 
the hypothesis of a dissipationless, merging formation 
for NSCs. We now wish to test whether the merger of 
globular clusters can result in a NSC similar to that of the 
MW. As a first step one can derive the globular cluster 
parameters required to give a peak density, pnsCj equal 
to the observed density of the NSC in the MW, and see 
whether such parameters are reasonable when compared 
to typical globular cluster properties Q- 

The stellar cluster central phase space density is given 
by 

PL -J. L_ na\ 

o K AnGo K r K ' U8J 

where pk is the cluster core density. After inspiral, we 
require Pnsc/ct nsc < Pk/<J K , or 
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so that 
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(20) 



The observed NSC density within 0.5 pc at the Galac - 
tic Center (GC) is p N sc ~ 10 6 M Q pc- 3 (jMerritt H20ifH ). 
Comparing this with Equation (|20p . we obtain 



ok 



lOkr 



"TK 

lpc 



< 15 



(21) 



This crude calculation shows that for a MW-like 
galaxy, the globular cluster parameters required to give 
the observed peak density of the NSC are quite reason- 
able. 

3.2. NSC formation via cluster migration 

A simple analytical model is developed in what fol- 
lows to evolve a population of stellar clusters subject to 
migration and dissolution in a galactic bulge and to cal- 
culate the influx of migrating clusters into the center of 
the galaxy. In this way we can address the possibility 
that a substantial fraction of the NSC mass in galaxies 
like the MW could have been assembled through cluster 
migration and mergers. 

1 I am indebted to D. Merritt for suggesting this calculation. 
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Fig. 1. — Formation of a NSC in a galactic bulge via cluster migration. Dashed and solid lines give respectively the density profiles of the 
NSC after 10 10 yr and the density of the background galaxy. Dot-dashed lines are the sum of galaxy and NSC density profiles. Red-dotted 
lines give t he power-law densit y model p(r) = 1.5 X 10 5 A/q (r/1 pc) - 2 , representative of the observed radial distribution of stars in the 
MW NSC (Schodcl ct al. 2009). We show results assuming that the total mass in stellar clusters is initially 10 per cent of the total galaxy 
mass and for M. = 4 X 10 6 M0(left panels), and M. = (right panels). The high mass truncation of the CIMF is M max = 10 7 Mq in the 
upper panels and M max = 10 6 Mq in the lower panels. Both large values of M. and small values of A-f max tend to reduce the excess density 
due to cluster infalls. This figure clearly demonstrates how MBHs can control the structure of NSCs forming through cluster migration 
and merging. 



We assume that the central properties of a globular 
cluster (i.e. ok and tk ) remain unchanged during inspi- 
ral and that r t > r^, where r t is the cluster tidal (limit- 
ing) radius given by (King 1119621) 



ok ( 3 d4> 
y/2 \r dr 



AirGp 



-1/2 



(22) 



with <fi the galactic potential, and a a "form factor" that 
depends on the density distribution within the cluster. 
For a King model with a large central concentration, 
awl. Equation (|22|) includes the tidal force due to 
the radial gradient in the galaxy potential, the centrifu- 
gal force from the cluster's orbit and assume that the 
density of the cluster goes to zero roughly at the radius 
where the force acting to remove a star is balanced by 
the attracting force from the rest of the cluster. 

The tidal radius of a King model orbiting in a galaxy 
with a power-law density profile, p(r) — po {r/r$) 7 , and 



containing a MBH at its center, is 
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The radius at which the tidal force from the MBH starts 
to dominate the tidal force from the stellar cusp (in 
galaxies containing both) is 



3(3 - 7 )M. 

747rp fo 

A King model also satisfies the relation: 



Gm t 



°K r t 
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(25) 



where m t is the truncated mass of the globular cluster 
whose radius is limited by the external tidal field. 

Given the cluster central velocity dispersion, Equa- 
tions ([23]) and ([25| can be combined to evaluate, at any 
radius, the cluster mass permitted by the galaxy tidal 
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field and the mass dispersed along the orbit. This cor- 
responds to a density enhancement with respect to the 
galactic background of 



Mr) 



K 



dn 



8irGr 2 dr 



(26) 



The resulting radial dependence of the density profile of 
stars in the growing NSC is easily found to be 



Ap(r) oc a 




cusp , 
black hole 



(27) 



Steeper density cusps in the distribution of background 
stars give lower values of the density slope for the result- 
ing NSC. Larger values of M, also give smaller densities 
near the center since the cluster models start being trun- 
cated at larger radii and their mass is dispersed over 
larger spatial scales. 

3.2.1. Dynamical Friction 
We obtain the cluster orbits by 



dard Chandrasekhar' s 
mula dChandrasekhar I J l943) 



using the stan- 
dynamical friction for- 
for a self-gra vitating 



cusp fe.g. JMerritt et al. 112001 lJust et al~| [2010): 
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The coefficient F is a weak function of 7, with F — 
(0.193, 0.302, 0.427) for 7 = (1, 1.5, 2). We set an 
initial limiting radius of 40 pc. If r t > 40 pc (i.e. the 
model is not truncated) , integrating Equation (|2"5|) yields 



r(t) = 



V (3- 7 ) 3/2 (6- 7 ) 
(4-7) 
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which gives a timescale to reach the center of the galaxy: 

(4-7) 



= 10 x 10 iu yr 
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with p ,5 = Po/5 M pc 3 , r ,7oo = r /700 pc, lnA 3 
In A '.') and m c \ t e = rn c \/ 1Q 6 M@. The coefficient depend- 
ing on 7 is approximately 1. 

After a cluster starts being truncated by the galactic 
tidal field, its mass also becomes a function of radius. 
In this case, we computed the cluster orbit by setting 
m c \ — m t in Equation f|28f>. obtaining 



r(t) 



(3- 7 ) 3 



(4- 7 )V2^4^Gp 



F(i) InAafxt 



possible presence of a MBH . Equation (l3~Tj) implies 
that the massive object comes to rest at the center of 
the stellar system in a time 



.(31) 



This equation takes into account mass loss due to the 
interaction with the galactic tidal field, but ignores the 



r* = 3 x 10 10 yr 



(4-7)^27- 
(3 -7)^(7) 



(32) 



In A 3 Vcsro^oo^Xio 



ro 
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with ax, 10 = c^/lO km s 1 and the coefficient depend- 
ing on 7 equals to (2.8, 4.27, 9.4) for 7 = (1, 1.5, 2). 

These basic formalism assumes that mass loss from the 
clusters is entirely due to their interaction with the exter- 
nal tidal field and it omits mass loss (evaporation) due to 
internal dynamics. Since the only clusters that can con- 
tribute to the growth of a NSC are very massive systems 
with relaxation times longer than their dynamical fric- 
tion timescale this simplification is quite reasonable and 
does not alter our results in any important way. How- 
ever, dynamical evaporation due to internal dynamics is 
also accounted for, in the sense that we do not follow the 
evolution of clusters which are disrupted (evaporated) on 
timescales shorter than the relevant dynamical friction 
timescales. 

3.2.2. Formation of the MW NSC 

Luminosity profiles of galaxies are well approximated 
by power laws with 1 < 7 < 2 at radii smaller 
than the effective radiu s of the stellar spheroid (R e ff, 
iTerzic & Graham II2005T ) . Since for MW-like galaxies the 
only clusters that can reach th e GC in one Hubble time 
are those initially at r < R e f / (|Milosavlievic 1120041) , and 
that within these radii the baryonic matter dominates the 
galactic potential, we represent the galaxy bulge by using 
a simple power-law density model: p(r) = po(r/ra) 1 , 
with po = (3 - 7)M snh /4 7rr 3 B M s „ h = 10 10 Afr, , 
r = 700 pc and 7 = 1 (e.g., McMilla nfc Dehnen |[2007^ . 
The scale length ro, is related to the bulge effective ra- 
dius via Reff/ro = (1.8,1.5,1) for 7 = (1,1.5,2). We 
assume that the clusters have i nitially the same distribu- 
tion o f stars in the galaxy (e.g. JAgarwal &: Milosavljevic I 
120111 ) and we assign their masses using the cluster ini- 
tial mass function (CIMF) , An/dM oc M~ 2 (|Bik et al. I 
120031 : Ide Grijs et al. 1 120031 ). and limiting mass values of 
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The vast majority of stars in a galaxy form in young 
protoclusters embedded within giant molecular clouds, 
but, due to expulsion of gas by feedback from su- 
pernovae and massive-star wind, only a small frac- 
tion ( ~ 5 %) of these systems survives the gas embedded 
phase (jTutukov Ill978l : lLada et al. 111984ft . The fraction 
of all star formation that occurs in stellar clusters, at 
the onset of the gas free-evolution, is usually referred to 
as "cluster formation efficiency". The cluster formatio n 
efficiency in the MW is - 10 % (|Lada fc Lada ll200l . 
Accordingly, in our model we assume that the total mass 
in clusters is initially M c \ = 0.1 x M spil . 

2 These equations also assume circular orbits. This is consis- 
tent with the well kno wn effect of orbital circularization due to 
dynamical friction (e.g., Cascrtano ct al. 1987; Hashimoto ct al. 
[2003 ). 

3 T his expression assumes that the density follows a IDehnenl 
(1993) profile at large radii, i.e., p(r) ~ r -4 for r 2> ro- 
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Fig. 2. — Effect of a galactic center MBH on the NSC properties. We assume M max = 10 7 Mq and that the total mass in clusters is 
initially 0.05 X M sp h (open circles), 0.1 X Af sp h (filled circles), or 0.3 X M sp h (star symbols). The properties of the galaxy are changed with 
M. according to the observed scaling relations between MBHs and the properties of their host spheroid. In the left panel we plot the ratio 
M./(M. + Moo) as a fun ction of M, , where M<3q is the mass deposited in the inner 30 pc of the galaxy. Dashed line is the observed 
correlation: Equation (1) in Graham & Spitlcr (2009). This plot measures the relative importance of M, to the total mass in the nuclear 
components for different black hole masses. The middle panel gives the importance (in terms of mass) o f the nuclear component ( MBH + 
NSC) relative to the host spheroid mass. Dashed line is the observed correlation given by Equation (2) in Graham & Spitlcr (2009). As an 
example, the right panel displays the density profile of background galaxy (solid lines), NSC (dashed lines) and the sum of them (dot-dashed 
lines) for an initial total mass in clusters equal to 0.1 X M B p& and for M, = 10 6 and 10 s Mq. Clearly, larger MBHs correspond to NSCs 
with lower densities. For M, > 10 s , the NSC central densities remain below the density of the galaxy at all radii. 



We start by identifying a "dissolution time" for glob- 
ular clusters due to dyna mical evaporation as a func- 
tion of their initial mass (IBaumgardt fc Makino 1 120031 : 
IGieles fc Baumgardt Il2008fc lLamers et al. 1120101 ): 

with P — 0.7, m c \ — rrit if the cluster is truncated, 
to = 0.3 x \dfl/d\n r| _1 and il the angular velocity 
at the galactocentric radius r. We then remove clus- 
ters with total lifetime, tdis/7, shorter than the dynam- 
ical friction time-scale, T*. This procedure reduces the 
initial total mass in globular clusters from 10 9 Mq to 
- 5 x 10 7 M & for M max = 10 6 Mq, and to - 2 x 10 8 Mq 
for M max = 10 7 Mq. By comparing Equation (|3"2")) with 
Equation (|3"3"|) one can easily show that if 7 < 3/2 and 
tdia/l > initially, this latter condition remains satis- 
fied during inspiral. 

We set the core radius of the globular clusters to be 
tk = 1 pc, roughly equal to the median value of the 
core radii listed in the Harris's compilation ((Harris II1996D 
of Galactic globular clusters. If a cluster reaches its 
tidal disruption radius its core mass is distributed uni- 
formly at the radius of disruption over a region of extent 
3 x rj{. When a cluster enters the inner 1 pc, we redis- 
tribute its mass using a Plummer sphere model of total 
mass mt(lpc) and core radius tr. To approximately ac- 
count for the mass that is progressively accumulated into 
the nucleus, the orbital decay of clusters was computed in 
order of increasing migration time. Then, for any inspi- 
ral, we added the quantity 3GM acc (< r)/r 3 to the terms 
in square brackets in Equation (|23|). with A/ acc (< r) the 
accumulated mass within r due to clusters with shorter 
orbital decay times. 

Finally, we computed the NSC density at any radius by 
summing up the contribution, Ap(r), of all clusters that 
within 10 10 yr reached that radius. After this time the 
total mass left in stellar clusters is ~ 10 7 Mq, almost in- 
dependent on the initial value of M max . The initial mass 



of the globular cluster system is therefore reduced from 
10 % to 0.1 % of the total Bulge mass after one Hub- 
ble time due to dynamical dissolution of the less massive 
systems, and also due to massive clusters that inspiral 
into center of the galaxy and dissolve locally to form the 
stellar nucleus. Figure Q] gives the results of such calcu- 
lations. 

In agreement with lAgarwal fc M ilosavli evic I (|2011[ ) we 
found that the mass of the forming NSC is mostly 
(> 90%) composed of clusters with initial masses > 
0.1 x M max , suggesting that the low mass clusters scarcely 
contribute to NSC formation in galaxies. Instead, our re- 
sults very much depend on the value of M max , since only 
very massive clusters can arrive in the central regions of 
the galaxy in a reasonable time without being destroyed 
by the galactic tidal forces in the process. 

In galaxy models without a MBH, NSC formation is 
more efficient, and even for M max = 10 6 Mq a NSC 
forms in the central few parsecs of the galaxy. If a MBH 
is present at the center of the galaxy, star clusters that 
come closer than r<ji S r from the center are disrupted and, 
as a result, the NSC density is limited inside this ra- 
dius (left panels in the figure). For a MBH mass of 
4 x 10 6 M Q and M max = 10 6 Mq, no clusters penetrate 
radii smaller than ~ 7 pc. From the lower-left panel 
in Figure Q] we can see that, in this case, the presence 
of a MBH will strongly inhibit the formation of a NSC. 
Taking M max = 10 7 M Q (upper-left panel) results in a 
smaller radius of the NSC core and higher central den- 
sities, since the most massive clusters can reach smaller 
galactocentric radii (~ 2 pc). 

A model that starts with a high mass truncation 
of M max = 10 7 M Q reproduces exceptionally well the 
observed mass density profile of stars in the Galactic 
NSC (red-dotted lines in the figure) outside 3 pc. We 
conclude that a scenario in which a large fraction of the 
mass of the Milky Way NSC is due to inf ailing globular 
clusters is in good agreement with current observational 
constraints. 
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The results presented here should be interpreted with 
caution when comparing the details of the NSC density 
profile to observations. Our modeling includes a series 
of simplifications that can have some impact on the in- 
ner structure of the resulting NSC. For example, consid- 
eration of the internal dynamics and mass spectrum in 
the globular clusters could enable mass segregation which 
can increase <7k and then allow the cluster stars to reach 
much smaller radii. A realistic treatment of these effects 
will require careful A-body simulations and is beyond the 
scope of this paper. Nevertheless, we note that the re- 
sults of our computations that inc lude a central MBH are 
very similar to those obtained in lAntonini et al. I (|2012[ ) 
via A-body simulations. In both cases, the density pro- 
file that results after the final inspiral event is charac- 
terized by a core of roughly the tidal disruption radius 
of the most massive clusters, ~ 3 pc, and an envelope 
with density that falls off as p ~ r~ 2 . These proper- 
ties are similar to those of the MW NSC, with the ex- 
ception of the core size, which in t he MW is somewhat 
smaller (~ 0.5 pc). iMerritt I (|2010() showed that such a 
core will shrink substantially via gravitational encounters 
in a time of 10 Gyr as the stellar distribution evo lves to- 
ward a Bahcall- Wolf cusp (jBahcall fc Wolf lfl976h . In our 
computations cluster inspiral occurs more or less contin- 
uously over the lifetime of the galaxy. The core resulting 
from the combined effects of cluster inspiral and relax- 
ation would therefore be somewhat smaller than that we 
found above, and closer to the observed core size (we 
discuss this point in more details below in H6.1j) . 

3.2.3. The role of MB Hs 

Figure [5] shows the results of additional computa- 
tions that explore the effect that varying Al, has 
on the structure of the NSC. In these plots we set 
M max = 10 7 Mq and we relate the mass of the 
galaxy to M. through the relation log [M,/M & ] ss 
8.4 + 1.91og [A/ S p h /7 x 10 10 M Q ] (IGraham I I2012bl) . 
The galaxy velocity disp ersion, a, is derived from 
the M. - a relati on (jFerrarese fc Merritt I 120001: 
IGebhardt et al. I I2001D . while the model scale length, 
r , is obtained from the effective radius: R e ff ~ 
1 kpc (M sph /10 10 M Q )/(o-/100 km s" 1 ) 2 , where here the 
normalization reproduces approxima tely the observed ef- 
fective radii of elli ptical galaxies (jForbes et al. 1 120081 : 
IGraham fc Worlev 1 120081 ) . In the left panel of Figure [2] 
we show the increasing dominance of the central MBH 
over the NSC stars for more massive galaxies. Our re- 
sults agree reasonably well with studies of galaxies at 
high mass end with M. > 10 8 Mq, and which ex- 
clude the presence of stellar nuclei in such systems (e.g., 
IHaring fc Rix 1 120041) . The middle panel of the figure 
displays how the total central mass in the nuclear com- 
ponent (i.e., MBH + NSC) divided by the stellar mass of 
the host spheroid varies with this latter quantity. Also 
in this case our results agree reasonably well with the 
observed correlation (dashed line in the plot). 

In order to reproduce the NSC masses obtained from 
observations, our model would require a high mass trun- 
cation of the CIMF of M max = 10 7 M Q and that < 10 % 
of stars in galaxies originate in stellar clus ters. Such a 
result is in agreement with the conclusion of lLeigh et al. I 
(|2012D that the globular cluster infall model strongly 



w o 

Q. — 



- 1 1 1 1 III 1 


■ i i 


1 I I I I I ■ | 


1 1 1 1 1 114 






y=0.5 
1 








1.5 








2 




f 
























1 1— 


1 L 


i 





10 

r(pc) 




n r 



■0 



6 7 8 

log M NSC (M Sun ) 



Fig. 3. — Upper panel gives the density profile of NSCs af- 
ter 10 10 yr in galaxy models with inner density profile slope 
7 = (0.5, 1, 1.5, 2). In these computations we set M sp ^ = 10 10 Mq, 
R e j / = 1 kpc, M, = 0, and M c i = 0.1 X M sph . The high mass trun- 
cation of the CIMF is M max = 10 7 Mq. The largest NSC central 
densities are found in models with small values of 7. In the middle 
and lower panels the properties of the NSC models are compared 
to those of real NSCs in early typ e galaxies (dots are data from 
ICote et al.1l200el; ISeth et al. I [2003) . The middle panel plots the 
nuclei in the size ratio versus mass ratio plane. Open circles corre- 
spond to the NSC models of the top panel with 7 = (0.5, 1, 1.5), 
where the relaxation time and effective radius of the NSC decrease 
with increasing 7. The blue-filled squares represent the model in 
the upper right panel of Figure [l] Lines indicate the critical value 
of R/R e ff over which nuclei expand for galaxy models described 
by Einasto indices n = 2 (continue line), n = 3 (dashed line) and 
n = 4 (dot-dashed line), (for comparison see Figure lc in IMerritt I 
120091 . The lower panel shows the dependence of NSC half-mass 
relaxation time on its total mass. 



under-predicts the observed NSC masses when assuming 
the present day number of globular clusters in galaxies. 
This assumption might be incorrect however, given that 
the number of globular clusters in a galaxy dr ops rapidly 
with time due to dynamical disruptions (e.g.. | Vesperinil 
H997L [19981 iFall fc Zhang1l200~Ll iGieles 1 120091 ) . In ad- 
dition, the cluster formation efficiency also varies sub- 
stantially with cosmic time, peaking at early epochs, in 
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gas-rich disk galaxies (pKruiisse n et al. I [2012). In fact, 
about 30 — 50 per cent of all stars in the Uni verse could 
have been formed in bound stellar clusters (|Kruiissen I 
120121) . 

The right panel of Figure [5] gives the density profile 
of a NSC that forms around a MBH of mass M, = 

10 6 — 10 8 Mq assuming an initial mass in clusters M c \ = 
0.1 x A/ sp h. Smaller black hole masses correspond to 
smaller tidal disruption radii for the infalling globu- 
lar clusters. As a consequence, the NSC peak density 
becomes progressively larger and the NSC core radius 
smaller as M, decreases. The mass delivered in the in- 
ner region of the galaxy is also a strong function of M.. 
For M. > 10 8 M Q , very little mass (< 10 6 M Q ) can be 
deposited in the inner ~ 30 pc of the galaxy. In this case 
the NSC density profile remains below the stellar density 
of the background galaxy. 

To summarize the results presented in this section, if 
we postulate that NSCs grow slowly through globular 
cluster migration and merging, then in order to form a 
NSC similar to that of the MW, our model will require 
M max ~ 10 r M Q and that ~ 10 per cent of the Bulge 
mass originated in stellar clusters. We showed that this 
fraction reduces to roughly 0.1 % of the total Bulge mass 
after 10 10 yr due to dynamical evaporation of the less 
massive clusters, and also because more massive clusters 
inspiral toward the galactic center and dissolve due to 
their interaction with the galactic tidal field. 

The presence of a pre-existing MBH at the center of the 
galaxy can have a strong impact on the inner structure of 
the forming nucleus. If a MBH of mass M, ~ 10 6 Mq sits 
at the center of the galaxy, massive globular clusters are 
disrupted at radii of order ~ 1 pc and the NSC density 
profile will have a core of roughly this radius. The cen- 
tral peak density of the resulting NSC progressively de- 
clines as M, increases. This suggests that NSCs forming 
around MBHs more massive than about 10 s Mq could 
have such low densities that they cannot be observed as 
distinct galactic components. 

4. DEPENDENCE ON 7 

In the above discussion we have considered Dehnen's 
models with inner density profile slope d\ogp/d\ogr = 
— 1 . In Figure [3] we relax this assumption and explore 
the formation of NSCs in galaxy models with different 
values of 7. In these computations we adopt Af max = 

10 7 Mq and M c \ = 0.1 x A/ sp h. The structural param- 
eters defining the galaxy model are M sp h = 10 10 Mq, 
R e ff — 1 kpc, and M. = 0. 

On the basis of Equation ([2l)|) we would expect the fi- 
nal (after 10 10 yr) densities in the inner ~ 100 pc of the 
forming NSC to be lower in galaxy models with steeper 
density profiles, provided that the dynamical friction 
timescale of massive stellar clusters remains shorter than 
the Hubble time. The results illustrated in the upper 
panel of Figure [3] agree with this basic prediction, sug- 
gesting that, independently on 7, the sinking timescale 
of massive clusters in these models is short enough that 
they can always reach the center of the galaxy and form 
a compact nucleus. 

The middle and bottom panels in Figure [3] compare the 
structural properties of our model NSCs to the properties 



of NSCs i n early- type galaxi es that were found to be nu- 
cleated in lCote et ai~l (j2006). For each of the NSC mod- 
els an approximation of its effective radius is obtained 
as the effective radius of the best fitting Sersic model to 
the NSC projected density profile within a galactocentric 
radius of 30 pc. Given R, the total NSC mass, Mnsc> is 
then obtained as the total mass within a radius twice the 
effective radius. A useful reference time is the relaxation 
time computed at R. Setting InA = 12, and ignoring 
the poss ible presence o f a MBH, the half-mass relaxation 
time is (jSpitzer II1987D : 



i rh = 1.75 x 10 5 



[r h (pc)f 2 N'/ \ 
(m/M )V2 



(34) 



where N is the total number of stars and to is the mass of 
a single star (to = Mq in Figure [3]). From the figure it is 
evident that our NSC models have structural properties 
that are similar to those of real nuclei. 

In the absence of a central MBH the dynamical evo- 
lution of a nuclear cluster is a competition between core 
collapse, which causes densities to increase, and heat in- 
put from the surrounding galaxy, which causes the clus- 
ter to expand and densities to decrease ([Kandrup 1119901 ; 
IQuinlan I [19961: IMerritt I [20091) . In a double-Plummer- 
law galaxy model, the maximum size of a NSC of mass 
MNSc/M S ph = 10~ 3 in order to resist expansion is 
R/Reff « 0.02 (|Quinlan|[l996[ ). This simple condi- 
tion, when compared to the middle panel of Figure [31 
implies that the core-collapse time of our models are al- 
ways shorter that the time for the nucleus to absorb en- 
ergy from the rest of the galaxy. However, when using 
more realistic Einasto profiles to describe the surround- 
ing galaxy our NSC models appear to be close to the 
transition region between systems that are in the prompt 
core collapse phase and systems for which the evolution 
is driven by heating from the galaxy. Our models will 
therefore either expand or contract depending on the 
"compactness" of the surrounding galaxy. 

5. ORBITAL DECAY IN THE CORE OF A GIANT 
ELLIPTICAL GALAXY 

Stellar spheroids more luminous than ~ 10 10 5 Lq 
often host high mass black holes (M, > 10 9 Mq) at 
their center, while they show no evidence for nucle- 
ation. Rather, the density profile of such "giant" ellip- 
tical galaxies is ob served to be flat inside the influenc e 
radius of the MBH (jKormendv Hl987HLauer et al.ll2002j) . 
Two factors have been invoked to explain the observed 
lack of NSCs in such systems: (i) NSCs and MBHs 
grow in competition from the same gas reservoir (|Escala I 
120071 ). If a MBH forms first this can prevent, thorough 
its feedback, a NSC from growing if the gas accretion rate 
is sm aller than the Eddington rate ( Navakshin et al."l 
2009). (ii) Massive black hole binaries forming dur- 
ing the last galaxy-galaxy major merger destroyed their 
host NSCs by ejecting stars fr om the inner galactic re- 
gions (|Bekki fc Graham |[2010h . 

Observations seem to point against both scenarios. 
Picture (i) seems to break down for NSC-dominated 
galaxies, given the extremely low accretion rates ob- 
served in bulgeless galaxies and that at least few 
of them also con t ain ce ntral MBHs. The model of 
iBekki fc Grahaml (|2010f) is disfavored due to the fact 
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Fig. 4. — Upper panel gives the orbital evolution of clusters with 
different values of the central velocity dispersion in the core of a 
M87-like galaxy. The dynamical friction timescale for clusters in 
such models is very long due to the lack of stars that move slower 
than the local circular velocity inside rj n fl. Even after 10 10 yr the 
cluster orbital radii have hardly changed from their initial values. 
The lower panel gives the density profile of the galaxy, assuming 
that the clusters did have enough time to migrate to the center: 
dotted, dashed and dot-dashed curves show the results of the accu- 
mulation of 500 clusters with ax = 30; 20 and 13 km s _1 respec- 
tively. Red lines correspond to computations without a MBH at 
the center of the galaxy. The galaxy background model is shown 
as a solid black line. Vertical marks correspond to the orbital ra- 
dius of the cluster in the model with an MBH after 10 10 yr and 
starting from r; n = 300 pc. Arrows give the tidal disruption radii 
of the clusters due to the MBH. At radii smaller than these our 
computations including a MBH are no longer valid. 



that what sets NSC disappearance does not seem to 
be galaxy morphology. Rather, there is a limit to 
the MBH/NSC mass ratio that fixes a sharp transi- 
tion ^h2HLS§l3^i^§_wlth_r^C_ajid MBH-dominated galax- 
ies (jNeumaver & Walcher Il2012f ). 

Of course, it is possible that nuclei are absent in bright 
galaxies because some mechanism prevents them from 
forming in the first place, or because they did not have 
time to reform after they were destroyed by the scouring 
effect of binary black holes. In what follows, we show that 
NSCs might in fact be difficult to form in such galaxies, 
due to long dynamical friction time-scales of star clusters 
and also the little mass that the clusters can deliver to 



the center due to the strong MBH tidal field. In order 
to d o so, w e repeat a similar analysis to that presented 
in § 13.2.21 but for a galaxy like M87, a giant elliptical 
with no evidence for NSC. We stress that the problems 
discussed here and in § 13.2.21 are substantially different. 
The influence radius of the MBH in a galaxy like M87 
extends much further out (a few hundreds of parsecs from 
the center) than the tidal disruption radius of a massive 
globular cluster (~ 10 pc), while for the MW we found 
r disr i; Hnfi- In the former case, the clusters will therefore 
move in a stellar background whose velocity distribution 
is directly influenced by the presence of a MBH, which 
will strongly affect the timescale for cluster inspiral as 
we now show. 

The distribution of field-star velocities has the follow- 
ing form near a MBH: 



/(«*) = 



r( 7 + i) i 



r( 7 . 



27^3/2, 



, 2 7 



2vi 



\ 7-3/2 



(35) 



where the normalizing constant corresponds to unit to- 
tal number. This expression gives the local distribution 
of velocities at a radius where the circular velocity is 
v c = (GM./r) 1 / 2 , assuming that the density of field stars 
follows p(r) — /5(r/f)~ 7 . The phase space density is zero 
for > v csc = 2 1/2 v c . 

As 7 — > 1/2 the velocity distribution (|3"5|) becomes pro- 
gressively narrower, and for 7 = 1/2 all stars have zero 
energy; in other words, the number of stars with < v c 
goes to zero as 7 approaches 1/2. In the case of a test 
particle moving in a circular orbit the standard Chan- 
drasekhar's formula will therefore predict zero frictional 
force. It turns out that in this situation the frictional 
force must be computed using a more general formula- 
tion that also includes the contribution from stars mov- 
ing faster than the test particle: 



f ~ f( v *< v ) 1 f(v*>v) 
J df ~ J df ~r J df 



(36) 



-4nG 2 m c ip(r)— x (hrA J dv^ f (v*)v* 



where v is the velocity of the infalling cluster. The sec- 
ond term in parenthesis represents the frictional force 
produced by stars moving faster than the massive par- 
ticle. ./V — body experiments verify the accuracy of this 
formula (jAntonini fc Merritt II2012T) . 

The fraction of stars that move slower than the local 
circular velocity can be computed as 



'v,<v c 



r( 7 + i) 



(37) 



dx x^-^^VT^ 



0Fr( 7 -l/2) J 1/2 

This function varies smoothly from when 7 = 0.5 to 
0.5 when 7 = 2. For the fast moving stars an (ad-hoc) 
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approximation is 



'v,>v c 



In 



; 0.1721 + 0.52807 - 0.1812 7 2 
-529.1892 exp(-14.45677) . 



(38) 



For 7 > 2 , I Vi ,>v c ~ 0.53, and when 7 ~ 2 the contri- 
bution from the slow stars is of order ~ InA larger than 

the frictional force due to the fast moving stars. 

Given these expressions, and setting L — y/GM,r in 
Equation ([3]) we have 



dr 
dt 



87rVGpf 7 n » t 7- i 5/2— v 

x [InA I Vi<Vc + J v .>v ] m c \r 5/2 7 



M. 



3/2 



(39) 



This expression can be then used to compute orbits 
for both point-like and extended objects (e.g., a stel- 
lar cluster that experiences mass loss during inspiral) 
that move in a stellar cusp near a MBH. We identi- 
fied our model with the center of a galaxy like M87 . 
We adopted M. = 3 x 10 9 M Q (jMacchetto et al. 1119971 : 
IGehbardt et"aDl20TTI). a core velocity dispersion a h = 
278kms- 1 (|Young et al. 1 119781: lLauer et all 11992ft and 
we used the relation a 2 — 4:nGph(rh/3) 2 with 77, = 
600 pc to obtain the core density: pt = 35 M Q pc~ 3 . 
Taking 7 = 0.5, the density at f = r in fl = 300 pc is 
p = 50A/ Q pc- 3 . 

For a test particle of mass m c i, integrating Equa- 
tion ([39j) yields 



r(t) = 



-3/2 _ 8^Gpf 7 ( 7 -3/2) 
M. 3/2 

x [InA I v ,<v c + ^v*>v c ] m d x t 



(40) 



-3/2 



for 7^3/2, and 

r(f) = r in exp -r- — [InA I Vt<Vc + I Vt > Vc \m cl x t 



Ml' 2 

(41) 

for 7 = 3/2. 

Equation (j40|) corresponds to a characteristic decay 
time for angular momentum loss (for 7 > 3/2): 



r. = 4 x 10 8 yr 



[InA I v .<v c + Iv t >vJ 



7-3/2 

, ,:i „> —1—3/2 -1 /f iTl \7-3/2 
xM ;9.5P 5 r 300 ™cl,6 J 



(42) 



where M.,9.5 = M./3 x 10 9 M , p 5 o = p/50 M Q pc" 3 , 
and f 30 o = r/300 pc. We note that, for 7 < 3/2, r, 
becomes negative and, formally, Equation (1391) gives an 
infinite decay time to the center. In this latter case the 
dynamical friction timescale can be re-defined as the time 
required to reach a radius half of its initial value: r. x 
(1 - I/2 7 - 3 / 2 ), with r. from Equation g2l . 

In the case of extended objects, we assume that the 
central properties of the globular cluster (i.e. ok and tk) 
remain unchanged during inspiral. The cluster's limiting 
radius is given by Equation (|23[) . which permits express- 
ing the satellite mass as a function of radius. Setting 



m c \ = m t in Equation (|39p and assuming that the galac- 
tic potential is dominated by the MBH, we obtain 



r(t) 



7-3 



2 3 / 2 (7 - 3)7rpr 7 



V3 

x[lnA Iv„<v c 



(43) 



*]°K X * 



and 



r. = 5 x 10 9 yr [lnA/v * <v - +/v * >vJ - (I/2 7 - 3 - l) 
3 — 7 ' 



'^•,9.5^50 ^300 a K,10 J 



r*. \ 7 — 3 

' ?.r). \ ' 



(44) 



We used Equations (|4"0)l and (f4"3")l to trace the orbital 
evolution of massive stellar clusters in the core of our 
galaxy model. We take a limiting value for the mass of 
mjc = 4 x 10 6 Mq, corresponding to the non-truncated 
model, when the cluster is far from the center. The sim- 
ulated orbits (upper panel of Figured]) demonstrate that 
the dynamical friction time-scale in the core of bright 
elliptical galaxies is extremely long. Also very mas- 
sive clusters (ctk ^ 30 km s _1 ) starting well within the 
galaxy core (r"j„ ~ 300 pc), do not reach the center after 
one Hubble time. 

In the lower panel of Figure HI we use Equations (|23|) 
and (f26|) to compute the density profile of a NSC result- 
ing from the accumulation of 500 equally massive glob- 
ular clusters in the center of our M87 galaxy model. At 
radii smaller than rdi sr (marked by black arrows in the 
plot) the clusters are fully disrupted by the interaction 
with the MBH and our integrations are no longer valid. 
A^-body simulations show that the density profile of stars 
will be very flat or even declining toward the MBH inside 
raisr, as the stars accumulate near the radius of d isrup- 
tion (|Kim fc Morris 1120031: iFuiii et al. 1120091 f2010h . 

Red lines in the lower panel of the figure give the den- 
sity profile of the nucleus when the mass of the MBH is 
set to zero. In this case, most of the cluster mass is de- 
livered in the inner ~ 100 pc and consequently the NSC 
is much more centrally concentrated. After 500 inspiral 
events, the NSC appears to be distinct from the galaxy 
background density profile. Apparently, NSCs assembled 
around MBHs have a spatially more diffuse configuration 
and lower densities than those forming in galaxies with 
no MBH. We conclude that the formation of NSCs in gi- 
ant ellipticals might be inhibited by the presence of their 
central MBH for two basic reasons: the long dynami- 
cal friction timescale of massive objects in the galaxy 
core; clusters are disrupted by the strong tidal field of 
the MBH producing a merger remnant with density pro- 
file that rises only modestly above that of the background 
galaxy. 

6. DISCUSSION 

6.1. Long term evolution of NSCs around MBHs 

In the previous sections we showed that NSCs resulting 
from the accumulation of globular clusters around MBHs 
have lower central densities and larger cores than those 
forming in galaxies without MBHs. Such a result is con- 
sistent with the results obtained using more s ophisticated 
A^— body simulations ()Antonini et al. ||2012[ ). 
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Our analysis however did not account for the effects 
of internal dynamics occurring in the NSC during and 
after its formation. In a preexisting NSC, the presence 
of a MBH woul d inhibit core collapse, c ausing instead the 
formation of a IBahcall fe Wolf] ((1979) cusp on the two- 
body relaxation timescale, f ollowed by a sl ow expansion 
as stars are tidally disrupted (Merritt 2009) . Whether or 
not a Bahcall-Wolf cusp will form depends on the initial 
core size of the nucleus and on its relaxation time. 

For nuclei similar to the MW NSC, relaxation times are 
too long to assume that they have reached such a colli- 
sionally relaxed state, but they are still sufficiently short 
that gravitational encounters would substantially affect 
their structure over the age of the galaxy. Two body 
gravitational interactions will cause the central density 
of the nucleus to increase and the core to shrink as 
the stellar distribution evolv es toward the Bahcall-Wolf 
form (e.g.. lPreto et al.~ll2004T ) . The time evolution of the 
core radi us can be approxima tely described by the ex- 
pression (jAntonini et al. 112012ft 

r-corc(^) = 1-57 pccxp [t/0.25t reIax ] , (45) 

where t rc iax is the relaxation time computed at the ra- 
dius of influence of the MBH. For the MW, i rclax ~ 
20 Gyr and M. 4 x 10 6 M Q . A black hole of 
this mass corresponds to a core of ~ 3 pc (Figure 1). 
From Equation (|45[) . we see that gravitational encoun- 
ters occurring during and after the formation of the 
NSC will reduce the core size to ~ 1 pc after 10 10 yr. 
Such a value is more consistent with the size of the 
core observed in the distribution of late-type (old) stars 
in the inner parsec of the MW (iBuchholz et al. 1 120091: 
Do et al. 112001 iBartko et al. II2OI0T lYusef-Zadeh et al. I 
20111: INadeen et al. I !20l3T Whether NSCs in other 
"power-law" galaxies will turn out to have stellar cores 
similar to that of the MW remains to be seen. 

6.2. The distribution of stellar remnants near MBHs 

NSCs containing MBHs are expected to be the 
main birthplace of gravitatio nal wave (GW ) sources for 
space- based interferometers ((Hughes 1 12003ft . These in- 
clude the capture of stellar-mass black holes (BHs) by 
MBH s, also called " extreme mass-ratio inspirals " (EM - 
RIs; IMerritt et all 120111: lAmaro-Seoane et al. I 12012ft . 
Understanding systems like the MW NSC and their ori- 
gin is therefore crucial for making predictions about the 
event rates for low frequency gravitational wave detec- 
tors. 

Most of the EMRI rate calculations reported in the lit- 
erature are derived under the assumption that the galac- 
tic nucleus had enough time to reach a state of mass seg- 
regation, w hich implies a high density of BHs near the 
center (e.g. iFreitae et al. 1120061: IS opman fc A lexander 
2006; [Alexander fc Hopman 1 12009).' This appears to 
be in conflict with observations which suggest an un- 
relaxed s tate f or t he distribution of stars at the GC. 
IMerritt I ((2010ft and lAntonini fc Merritt I (|2012ft demon- 
strated that in the absence of an initial cusp in the stars, 
even after 5 — 10 Gyr the density of BHs could remain 
substantially below the densities inferred from steady- 
state models. 

All these previous studies, however, assumed that the 
stellar BHs had the same phase-space distribution ini- 



tially as the stars. This is most likely to be a poor as- 
sumption if a substantial fraction of the NSC mass comes 
from orbitally decayed globular clusters. In the merger 
model the resulting distribution of stellar remnants will 
reflect their distribution in their parent clusters just be- 
fore they reach the center of the galaxy. 

In a dense stellar cluster, BHs formed by the super- 
nova explosions of the most massive stars tend to segre- 
gate into the cluster core and form a sub-cluster of BHs 
whi ch dynamicall y decouples from the rest of the clus- 
ter ((Spitzer 111987ft . When the central density of BHs be- 
comes large enough, BH-BH binary formation becomes 
efficient. Subsequent dynamical interactions involving 
binary BHs and higher multiplicity systems will tend to 
eject the BH s from the cluster until on l y a few of them 
are le ft (e.g.. iDowing et al. Il2010l 120111: iBaiieriee et al. I 
120101) . 

We may consider two possibilities 0: (i) the dynamical 
friction timescale of star clusters is short compared to 
the encounter driven evaporation timescale of their BH 
sub-cluster. In this case, the mass-segregted BH sub- 
clusters, due to their high central den sities, can reac h 
galactocentric radii as small as ~ 0.1 pc ((Antonini 12 012). 
The preferential removal of stars from the outer parts of 
the clusters by the strong galactic tidal field might lead to 
the formation of a NSC with a central over-abundance of 
BHs when compared to predict i ons fr om standard mass 
functions (Ba neriee fc Kroupa 1 1 2 1 If) . (ii) The orbital 
decay timescale is shorter than the evaporation timescale 
of the cluster BH sub-system. Most of the BHs will be 
dynamically ejected in this case before the cluster reaches 
the center of the galaxy. As a consequence of this, very 
few BHs are delivered to the vicinity of the MBH. 

All of the above mentioned topics require a dedicated 
study which we defer to a future paper. 

6.3. In- situ star formation 

Our work shows that a purely dissipationless merger 
scenario can explain, without obvious difficulties, the ba- 
sic physical properties of NSCs. On the other hand, we 
note that star forming events triggered by rapid inflow of 
gas in the center of a galaxy may also contribute to the 
growth of its NSC. The fact that NSCs in general tend to 
have a wide range of stellar ages, including young stellar 
populations, seems to confirm the hypothesis that the 
nuclei formation hi stories were in par t governed by loca l 
or internal factors ((Rossa et al. Il2006t iSeth et al. ll2006ft . 

As an example, our GC contains a large population of 
young massive stars that most likely formed locally fol- 
lowing the infall and fragmentat ion of a dense gaseou s 
clump in the vicinity of Sgr A* (|Paumard et al. Il2006l) . 
Such bursts may occur continuously over the age of the 
Galaxy and thus produce a significant fraction of its NSC 
mass. Stellar population synthesis studies also show that 
the GC appears to have undergone continuous and re- 
current star formation over the last 10 Gyr, but it is not 
possible to fit the observations with ancient burst mod- 
els, such as would be appropriate for an ol d population 
of sta rs that originated in globular clusters ((Figer et al. I 
12001 . 

4 We make here the reasonable assumption that the timescale 
for a clusters to reach the center of the galaxy is long compared to 
the mass-segregation timescale of their BH population. 
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Previous work incorrectly used the above arguments 
to argue again st a dissipationless origin fo r NSCs (e.g., 
lMilosavlievi6ll2004tlNa"vakshin et al. 112009ft . It is impor- 
tant to stress that although observations probe recent 
and episodic star formation they do not exclude that the 
bulk of the GC stellar population is in old stars. In 
fact, the GC luminosity function appears to be consis- 
tent with a star formation history in which a large frac- 
tion (about 1/2) of the mass consists of old (~ 10 Gyr) 
stars and the remainder is d ue to continuou s star forma- 
tion (lAntonini et al. 1 12012). Accordingly. iPfuhl et al. 1 
(|2011ft found that about 80 per cent of the stellar mass 
in the inner parsec of the Galaxy is in old stars that 
formed more than 5 Gyr ago. We add that, although 
globular clust ers in our Galaxy are ex clusively old stellar 
systems (e.e.. iRosenberg et al. Ill999l ). this might not be 
always the case in other galaxies. It is therefore not clear 
what fraction of the mass in NSCs is contributed by local 
gaseous fragmentation. More investigation needs to be 
done in order to explore the implications of in-situ star 
formation as the origin of galactic nuclei. 

6.4. NSCs morphology and kinematics 

The morphology and kinematics of NSCs are of great 
importance for understanding their origin. 

Aspher ical NSCs a r e com monly observed in external 
galaxies. iSeth et al~1 ((2006) found that the three edge- 
on late-type galaxies IC 5052, NGC 4206 and NGC 4244 
have nuclei that are strongly elongated along the plane 
of their host galaxies disks. Such clusters show evidence 
for multiple morphological components, with a disk-like 
young stellar population superimposed on an older more 
spherical component. The radial velocity map of the 
nucleus in NGC 4244, the nearest of these three galax- 
ies (D = 4.1 Mpc), shows evidence for strong rotation, 
30 km s _1 at 10 pc from the center, co mpared to a 
centr al velocity dispersion of ~ 28 km s _1 (|Seth et al. I 
12008ft . There is also evidence for flattening and rotation 
in th e M33 nucleus (jLauer et al. 11199 8; Matt hews et al. I 
1999). The M33 NSC is elongated along the major-axis 
of the galaxy and rotates at ~ 8 km s" 1 , while the cen- 
tral v elocity dispersion is ~ 27 km s _1 (|Gebhardt et al. I 
12001ft . The MW NSC also appears t o be rotating paral- 
lel to the overall Ga lactic rotation (|Trippe et al. 1 120081 : 
iSchodel et aT~l 12009ft . Unfortunately, as a consequence 
of the strong interstellar extinction along the GC line of 
sight, our knowledge of the Galactic NSC morphology 
and size remains very limited. 

Is the observed rotation of NSCs consistent with the 
predictions of a dissipationless model for NSC forma- 
tion? Of course if a large fraction of the NSC mass 
was formed via accretion of clusters isotropically dis- 
tributed throughout the galaxy, no net rotation would 
be expected, since there will not be any preferred di- 
rection for inspiral. If the primary formation process is 
instead gas accretion from the galactic disk, this would 
naturally explain both the flattening and the fact that 
NSCs rotate parallel to their host galaxy rotation. 

However, we believe that the evidence for NSC rota- 
tion, at least in late-type galaxies, is not a strong ar- 
gument against the globular cluster merger model. For 
instance, clusters falling into the GC could have origi- 
nated in the inner part of the Galactic disk and they will 



therefore share in its rotation (Hartma nn et al.l 120 lift . 
Another possibility is that globular clusters crossing the 
Galactic disk experience a greater fri ctional force from 
the increased local stellar/gas density (jBekki 112010ft and 
hence they can be dragged down into the disk plane and 
transported into the central region of the galaxy where 
they then accumulate to form a dense nucleus. In ei- 
ther case, the forming NSC will appear to rotate in the 
same sense of the Galaxy. These hypotheses lead to basic 
predictions that might be testable with future observa- 
tions; specifically: (i) NSCs in early-type galaxies, which 
do not have extended massive disk-like structures, might 
have slower rotation with respect to NSCs of spiral galax- 
ies, (ii) There should be some mild depletion of stellar 
clusters in the Galactic disk as compared to off the disk 
plane 0. 

In conclusion, comparing kinematic data with simula- 
tions for distinguishing between gas and cluster accretion 
may be difficult, since the detailed structure of simulated 
NSCs varies in an important way with the orbital initial 
configuration of the infalling clusters. Future observa- 
tional work should be able to provide more reliable mod- 
els for the spatial distribution of globular clusters near 
the center of galaxies, where timescales for infall are rea- 
sonably short. 

6.5. The lack of NSCs in faint galaxies 

Galaxies fainter t han Mb ~ —12 do not contain promi- 
nent stellar nuclei (Ivan den Bergh 111986ft . The lack of 
NSCs in dwarf spheroidals was explaine d in the con- 
text o f the merger formation scenario in iGoerdt et al. I 
( 2006) . These authors used high resolution TV-body sim- 
ulations to study the orbital decay of globular clusters 
in dark matter haloes with a core-density structure in 
the ce ntral regions, similar to wh at inferred by observa- 
tions (|Prvor fc Kormendv 1 11990ft . In these models the 
standard Cha ndrasekhar's dynamical friction formula 
breaks down (Tremaine & Weinberg 1984; Weinberg 
1986) and one finds that a massive obj ect stalls at 
roughly the dark matter core radius (e.g., iRead et al. I 
[20061 llnoue 1 12009ft . This suggests that the lack of NSCs 
in the fainter spheroidals could be simply a consequence 
of the arbitrarily long sinking timescales of stellar clus- 
ters in these systems. 

7. CONCLUSIONS 

In this paper we considered a dissipationless formation 
model for NSCs where globular clusters orbital decay and 
merge at the center of a galaxy to form a compact nu- 
cleus. Our main results are summarized below. 

1 The observed scaling relation between NSC masses 
and the velocity dispersion of their host spheroids, 
the A/nsc — °~ relation, is difficult to reconcile with 
a purely dissipative formation model for the nu- 
clei. These models predict that Mnsc is a steeply 
rising function of a. The observed Mnsc — & rela- 
tion is instead in agreement with the predictions of 
a dissipationless formation model. Dissipationless 
formation modes produce relations that are sub- 
stantially shallower than the corresponding MBH 

5 I am thankful to A. Madigan for pointing this out. 
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scaling relations and are therefore more consistent 
with observations. 



2 The globular cluster merger model, in the absence 
of a central MBH, naturally reproduces the ob- 
served relation between the size of galactic nuclei 
and their total luminosity, R <x y/Cwsc- When a 
MBH is present, the dependence of the NSC ra- 
dius on its mass is substantially weaker than the 
observed relation because the size of the NSC is 
mainly determined by the fixed tidal field of the 
MBH. 



3 We derived explicit expressions for the orbits of 
globular clusters owing to dynamical friction and 
subject to mass-loss due to their tidal interaction 
with the galaxy (Equation I31[) or with a central 
MBH (Equation 1^3]) . These expressions were used 
to (i) address the possibility that NSCs in galaxies 
similar to the MW could have been assembled via 
cluster migration and mergers; and (ii) to follow the 
orbital evolution of globular clusters in the core of 
bright elliptical galaxies. 



4 NSCs that form through the mergers of globular 
clusters have a density profile characterized by a 
parsec-scale core and an envelope that falls off as 
p ~ r~ 2 . These properties are similar to those of 
the MW NSC. A NSC with mass comparable to 
that of the MW NSC is obtained by assuming an 
initial total mass in stellar clusters which is consis- 
tent with the cluster formation efficiency inferred 
from observational studies of embedded clusters in 
Galactic molecular clouds (Figure [1]). 



5 A pre-existing MBH at the center of the stellar 
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spheroid has a strong impact on the structure of the 
growing stellar nucleus (Figure [2]) . Tidal stresses 
from a MBH disrupt the clusters when they pass 
within their tidal disruption radius, rdi sr , limiting 
the density within that radius. Hence, the density 
profile of the resulting NSC will also have core size 
of ~ ^disr- Removing the MBH from the galaxy 
allows the stellar clusters to keep their inner struc- 
ture almost unchanged leading to the formation of 
a NSC with higher peak densities. We find that 
observing NSCs as a distinct component could be 
more difficult in stellar spheroids with more mas- 
sive MBHs. 



6 For globular clusters orbiting in the core of a mas- 
sive elliptical galaxy like M87, the timescale to 
reach the center is much longer than one Hubble 
time (Figure 0}. The essential reason for this is 
that the phase-space density of a shallow density 
cusp of stars around a MBH falls to zero at low en- 
ergies. In addition, when a MBH is included, the 
tidal field near n n fl becomes much stronger, result- 
ing in faster mass loss and fewer stars deposited to 
the center. Based on these facts, we conclude that 
the formation of compact nuclei in bright galaxies 
might be inhibited by the presence of their central 
MBHs. 
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